Quantum Monte Carlo Calculations with Chiral Effective Field Theory Interactions 
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We present the first Quantum Monte Carlo (QMC) calculations with chiral effective field theory 
(EFT) interactions. To achieve this, we remove all sources of nonlocality, which hamper the inclusion 
in QMC, in nuclear forces to next-to-next-to- leading order (N 2 LO). We perform Auxiliary-Field 
Diffusion Monte Carlo (AFDMC) calculations for the neutron matter energy up to saturation density 
based on local leading-order, next-to-leading order, and N 2 LO nucleon-nucleon interactions. Our 
results exhibit a systematic order-by-order convergence in chiral EFT and provide nonperturbative 
benchmarks with theoretical uncertainties. For the softer interactions, perturbative calculations are 
in excellent agreement with the AFDMC results. This work paves the way for QMC calculations with 
systematic chiral EFT interactions for nuclei and nuclear matter, for testing the perturbativeness 
of different orders, and allows for matching to lattice QCD results by varying the pion mass. 

PACS numbers: 21.60.Ka, 21.65.Cd, 21.30.-x, 26.60.-c, 02.70.Ss 
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Chiral EFT has revolutionized the theory of nuclear 
forces by providing a systematic expansion for strong 
interactions at low energies based on the symmetries 
of quantum chromodynamics [IHll- Chiral interactions 
have been successfully employed in calculations of the 
structure and reactions of light nuclei lMA, medium- 
mass nuclei and nucleonic matter 14-201. While 

continuum QMC methods are very precise for strongly in- 
~ 2 



including neutron matter [25 



uation using propagators of the form: 
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teracting systems 

[29l |. and have provided pioneering calculations of light 
nuclei 
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24], QMC has not been used with chiral EFT 
interactions due to nonlocalities in their present imple- 
mentation in momentum space. In this Letter, we take up 
this challenge and combine the accuracy of QMC meth- 
ods with the systematic chiral EFT expansion. As an ap- 
plication, we study the neutron matter equation of state 
at nuclear densities. Neutron matter constitutes an excit- 
ing system because of its connections to ultracold atoms 
and its importance for the physics of neutron-rich nuclei, 
neutron stars, and supernovae. 

First, we explain how to remove all sources of nonlo- 
cality in chiral EFT interactions to N 2 LO and present 
local nucleon-nucleon (NN) interactions at leading-order 
(LO), next-to-leading order (NLO), and N 2 LO based on 
Ref. [30(. We use the developed chiral potentials for the 
first time in QMC calculations to study neutron matter 
order-by-order including theoretical uncertainties. The 
nonperturbative QMC results provide many-body bench- 
marks and enable us to test perturbative calculations for 
the same interactions. 

The difficulty of handling nonlocal interactions in 
QMC methods results from how interactions enter. Con- 
tinuum QMC methods are based on a path-integral eval- 



G(R,R') = (R| 



-StO 



R'), 



(1) 



where R = (ri, r% ■ . . r/y) is the configuration vector of 
all N particles (plus spins and other quantum numbers) , 
St is a step in the imaginary-time evolution, and the op- 
erator O takes into account the kinetic energy and the 
interaction part of the Hamiltonian. The implementa- 
tion of continuum QMC methods relies on being able 
to separate all momentum dependences as a quadratic 
X)i=i Pi term, which is the case for local interactions, 
but not for general momentum-dependent, nonlocal in- 
teractions (spin-orbit interactions, linear in momentum, 
are manageable) . In the local case, the propagator for the 
momentum-dependent part is a Gaussian integral that 
can be evaluated analytically, and the effects of interac- 
tions only concern the positions of the particles. 

Chiral EFT interactions are based on a momentum 
expansion and are therefore naturally formulated in mo- 
mentum space [3, 0| • To regularize interactions at high 
momenta, one employs regulator functions, usually of the 
form f(p) = e~( p / A ) " and f(p'), where p = (pi — P2)/2 
and p' = (p^ — P2V2 are the incoming and outgoing 
relative momenta, respectively. Upon Fourier transfor- 
mation, this leads to nonlocal interactions V(r,r') al- 
ready due to the choice of regulator functions. The other 
sources of nonlocality in chiral EFT are due to contact 
interactions that depend on the momentum transfer in 
the exchange channel k = (p' + p)/2 and due to k- 
dependent parts in pion-exchange contributions beyond 
N 2 LO. In contrast, dependences on the momentum trans- 
fer q = p' p are local, and lead to nonlocalities only 
because of the regulator functions used. 
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To avoid regulator-generated nonlocalities for the long- 
range pion-exchange parts of chiral EFT interactions, we 
use the local coordinate-space expressions for the LO 
one-pion-exchange, NLO and N 2 LO two-pion-exchange 



interactions 



3l[ |32[ and regulate them directly in coor- 
dinate space using the function fi on g{r) = 1 — e~( r ' R °} , 
which smoothly cuts off interactions at short distances 
r < Rq, while leaving the long-range parts unchanged. 
So, Rq takes over the role of the cutoff A in momentum 
space. To regularize the pion loop integrals of the two- 
pion-exchange contributions, we^use a spectral-function 
regularization [12] with cutoff A = 800 MeV. For the 
N 2 LO two-pion-exchange interactions we take the low- 
energy constants c\ = — 0.81 GcV -1 , c 3 = — 3.4GeV~ , 
and C4 = 3.4 GeV -1 as in the momentum-space N 2 LO 



potential of Ref. [33 1. 

To remove the k-dependent contact interactions to 
N 2 LO, we make use of a freedom to choose a basis of 
short-range operators in chiral EFT interactions (simi- 
lar to Fierz ambiguities). At LO, one usually consid- 
ers the two momentum-independent contact interactions 
Cs + Ct <t\ ■ er 2 . However, it is equivalent to choose 
any two of the four operators 1, a\ ■ er 2 , T\ • T2, and 
<7i • er 2 Ti • T2, with spin and isospin operators tTj,Tj, be- 
cause there are only two S-wave channels due to the Pauli 
principle. It is a convention in present chiral EFT inter- 
actions to neglect the isospin dependence, which is then 
generated from the exchange terms. 

We use this freedom to keep at NLO (order Q 2 ) an 
isospin-dependent q 2 contact interaction and an isospin- 
dependent {cri ■ q)(<r 2 • q) tensor part in favor of a nonlo- 
cal k 2 contact interaction and a nonlocal (<xi ■ k)(er 2 • k) 
tensor part. This leads to the following seven linearly 
independent contact interactions at NLO that are local, 



V. 



NLO _ r 2 , n 2 _ _ 
short +<^2q Tl • T2 



C 4 q 2 Tl ■ T 2 ) (Ti ■ <T 2 
(cri + er 2 ) • q x k 



{C 3 q 

.c 5 

I — 
2 

C 6 {(Ti ■ q)(<T 2 • q) 

C 7 (o-i • q)(cr 2 • q) n • r 2 , 



(2) 



where the only k-dependent contact interaction (C5) is a 
spin-orbit potential. Because at NLO the only two pos- 
sible momentum operators allowed by symmetries are q 2 
and k 2 (or equivalently p 2 + p' 2 and p • p'), and simi- 
larly for the tensor parts, it is thus possible to remove all 
sources of nonlocality in chiral EFT to N 2 LO. In addi- 
tion, the leading 3N forces at N 2 LO can be constructed as 
local interactions [13, 111] , but we will first focus on QMG 
calculations with chiral NN interactions. The next-higher 
order (Q 4 ) NN contact interactions enter at N 3 LO, and 
there are too many possible operators involving k, so that 
they cannot be traded for isospin dependence completely. 
Therefore, chiral EFT interactions will contain nonlocal 



TABLE I. Short-range couplings for R = 1.2 fm at LO, NLO, 
and N 2 LO (with a spectral-function cutoff A = 800 MeV) [11. 
The couplings Ci_7 are given in fm 4 while the rest are in fm 2 . 
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terms at N 3 LO, but one may expect that these high-order 
nonlocal parts can be treated perturbatively. 

Upon Fourier transformation, these LO and NLO con- 
tact interactions lead to local smeared-out delta func- 
tions 5r (r) and their derivatives, when a local regulator 
/iocai(<Z 2 ) is used. We implicitly define the local regu- 
lator by taking <5R (r) ~ e-( r / R °) with an exponential 
regulator (with the same scale Rq) similarly as for the 
long-range parts. We thus have for the LO contact inter- 
actions in coordinate space 



1 C S ,T /localteV*' = C S . 



(27T) 



(3) 



where the prefactor is determined by normalization. The 
analogous local expressions involving the NLO contact in- 
teractions are obtained by replacing Cs.t with the seven 
different operators of Eq. (j2J|. Finally, for the range of 
the scale Rq we consider Rq = 0.8 — 1.2 fm corresponding 
to typical momentum cutoffs A ~ 600 — 400 MeV in chiral 
EFT interactions. 

The low-energy couplings Cs,t at LO plus Ci_7 at 
NLO and N 2 LO are fit in Ref. ' 30] for different Rq to 
the NN phase shifts of the Nijmcgen partial-wave analy- 
sis [3^] at laboratory energies E\ a \, = 1, 5, 10, 25, 50, and 
100 MeV. The reproduction of the isospin T = 1 S- 
and P-waves is shown order- by-order in Fig. [I] where the 
bands are obtained by varying Rq between 0.8 — 1.2 fm 
and provide a measure of the theoretical uncertainty. For 
the R = 1.2 fm N 2 LO NN potential, we list the low- 
energy couplings at LO, NLO, and N 2 LO in Table H At 
N 2 LO, an isospin-symmetry-breaking contact interaction 
{Cnn for neutrons) is added in the spin 5 = channel (to 
Cs — 3Cr), which is fit to a scattering length of —18.8 fm. 
As shown in Fig. [TJ the comparison with NN phase shifts 
is very good for Ej a b < 150 MeV. This is also the case for 
higher partial waves and isospin T = channels, which 
will be reported in a following paper. In cases where 



3 




Lab. Energy [MeV] Lab. Energy [MeV] Lab. Energy [MeV] Lab. Energy [MeV] 

FIG. 1. (Color online) Neutron-proton phase shifts as a function of laboratory energy Ei^b = 2p 2 /m in the 1 So, 3 Po, 3 Pi, and 
3 P2 partial waves (from left to right) in comparison to the Nijmegen partial-wave analysis (PWA) [3(|. The LO, NLO, and 
N 2 LO bands are obtained by varying Ro between 0.8 — 1.2 fm (with a spectral-function cutoff A = 800 MeV). 



there are deviations for higher energies (such as in the 
3 P 2 channel of Fig. [T]), the width of the band signals sig- 
nificant theoretical uncertainties due to the chiral EFT 
truncation at N 2 LO. The NLO and N 2 LO bands nicely 
overlap (as shown for the cases in Fig. [1]), or are very 
close, but it is also apparent that the N 2 LO bands are 
of a similar size as at NLO. This is because the width of 
the bands both at NLO and N 2 LO shows effects of the 
neglected order-Q 4 contact interactions. 

We then apply the developed local LO, NLO, and 
N 2 LO chiral EFT interactions in systematic QMC cal- 
culations for the first time. Since nuclear forces con- 
tain quadratic spin, isospin, and tensor operators (of the 
form <t\ A 1 ^ cr?) i the many-body wave function cannot 
be expressed as a product of single-particle spin-isospin 
states. All possible spin-isospin nucleon-pair states need 
to be explicitly accounted for, leading to an exponential 
increase in the number of possible states. As a result, 
Green's Function Monte Carlo (GFMC) calculations are 
presently limited to 12 nucleons and 16 neutrons [jij ]. In 
this Letter, we would like to simulate 0(100) neutrons 
to access the thermod yna mic limit. We therefore turn to 
the AFDMC method [37|, which is capable of efficiently 
handling spin-dependent Hamiltonians. 

Schematically, AFDMC rewrites the Green's function 
by applying a Hubbard-Stratonovich transformation us- 
ing auxiliary fields to change the quadratic spin-isospin 
operator dependences to linear. As a result, when ap- 
plied to a wave function that is a product of single- 
particle spin-isospin states, the new propagator indepen- 
dently rotates the spin of every single nucleon. Using this 
approach, central and tensor interactions can be fully in- 
cluded in an AFDMC stochastic simulation. For the case 
of neutrons, it has also been possible to include spin-orbit 
interactions in AFDMC 



We first studied finite-size effects and the dependence 
on the Jastrow correlations in the trial Jastrow-Slatcr 
wave function (in continuum QMC there are no dis- 
cretization effects). The dependence on particle num- 
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FIG. 2. (Color online) Neutron matter energy per particle 
E/N as a function of density n calculated using AFDMC 
with chiral EFT NN interactions at LO, NLO, and N 2 LO. 
The statistical errors are smaller than the points shown. The 
lines give the range of the energy band obtained by varying 
Ro between 0.8 — 1.2 fm (as for the phase shifts in Fig. [1]), 
which provides an estimate of the theoretical uncertainty at 
each order. The N 2 LO band is comparable to the one at NLO 
due to the large a couplings in the N 2 LO two-pion exchange. 



ber was found to be nearly identical to that of the non- 
interacting Fermi system, consistent with results using 
phenomenological NN potentials [39(. Therefore, we per- 
formed calculations for an optimal number of 66 parti- 
cles, while also including contributions from the 26 cells 
neighboring the primary simulation box. We also com- 
pared the neutron matter energy at a density 0.1 fm -3 
starting from no to full Jastrow correlations based on the 
same Rq local chiral NN interactions versus Jastrow cor- 
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FIG. 3. (Color online) The AFDMC N 2 LO band of Fig. g] 
in comparison to perturbative calculations of the neutron 
matter energy using the same local N 2 LO NN interactions. 
The lower (upper) limit of the AFDMC N 2 LO band is for 
Ro = 1.2 fm (Ro = 0.8 fm), corresponding to a momentum 
cutoff A ~ 400 MeV (A ~ 600 MeV). Perturbative results are 
shown for Hartree-Fock plus second-order contributions (2nd 
order) and including third-order particle-particle and hole- 
hole corrections (3rd order). The bands at 2nd and 3rd order 
are obtained by using a Hartree-Fock or free single-particle 
spectrum. For the softer Ro = 1.2 fm interaction (narrow 
purple bands), the third-order corrections are small and the 
perturbative third-order energy is in excellent agreement with 
the AFDMC results, while for the harder Ro = 0.8 fm interac- 
tion (light red bands), the convergence is clearly slow. At low 
densities, we also show the QMC (2010) results of Ref. [29| . |43T |. 



relations of the hard Argonne v' a potential, as a first step 
in probing the general dependence on the Jastrow term. 
For the softer Ro = 1.2 fm (Rq — 0.8 fm) interactions the 
changes of the energy per particle are at most 0.1 MeV 
(0.6 MeV), which corresponds to 1% (5%) changes. This 
appears to be related to the way the propagator is sam- 
pled with tensor and spin-orbit interactions and will be 
studied in detail in a forthcoming paper. The exact re- 
sults should be independent of the trial wave function, 
but we consider Jastrow correlations based on the same 
Rq interactions more consistent and use these. 

In Fig. [5] we present first AFDMC calculations for the 
neutron matter energy with chiral EFT NN interactions 
at LO, NLO, and N 2 LO. Our results represent nonper- 
turbative energies for neutron matter based on chiral 
EFT beyond low densities. For neutrons, the AFDMC 
method has been carefully benchmarked with nuclear 
GFMC [Z(J, |4lj , with the fermion sign problem after re- 
leasing the fixed node having minimal effects. At each 



order, the full interaction is used both in the propagator 
and when evaluating observables. The lines in Fig.[2]give 
the range of the energy obtained by varying R between 
0.8 — 1.2 fm, where the softer R Q = 1.2 fm interactions 
yield the lower energies. At low densities (low Fermi mo- 
menta), as expected the energy is well constrained at LO, 
with small corrections at NLO due to effective range ef- 
fects [H, HU ■ AFDMC enables us to present results up 
to saturation density (and higher, but we emphasize that 
the contributions of 3N forces will become important for 
densities n > 0.05 fm -3 [16]). At LO, the energy has a 
large uncertainty. The overlap of the bands at different 
orders in Fig. [2] is excellent. In addition, the result that 
the NLO and N 2 LO bands are comparable is expected 
due to the large Cj entering at N 2 LO; this is similar to 
the phase shift bands in Fig.Q] Therefore, our first QMC 
results for neutron matter exhibit a systematic order-by- 
order convergence in chiral EFT. 

Our AFDMC results provide first nonperturbative 
benchmarks for chiral EFT interactions at nuclear densi- 
ties. We have performed perturbative calculations follow- 
ing Refs. [H, G3, El based on the same local N 2 LO NN 
interactions at the Hartree-Fock level plus second-order 
contributions and including third-order particle-particle 
and hole-hole corrections. At each order, we give bands 
obtained by using a Hartree-Fock or free single-particle 
spectrum. The perturbative energies are compared in 
Fig. to the AFDMC N 2 LO results. For the softer 
Ro = 1.2 fm (A ~ 400 MeV) interaction, the third-order 
corrections are small and the perturbative third-order 
energy is in excellent agreement with the AFDMC re- 
sults, while for the harder R = 0.8 fm (A ~ 600 MeV) 
interaction, the convergence is clearly slow. This is the 
first nonperturbative validation for neutron matter of the 
possible perturbativeness of low-cutoff A ~ 400 MeV in- 
teractions [Hj]. Finally, in the low-density regime, the 
results in Fig. [3] match on to the QMC calculations of 
Ref. ET 



431 ] based on central interactions that reproduce 
the large neutron-neutron scattering length and the ef- 
fective range physics. 

In summary, we have presented the first QMC calcu- 
lations with chiral EFT interactions. This was achieved 
by using a freedom in chiral EFT to remove all sources 
of nonlocality to N 2 LO. We have constructed local LO, 
NLO, and N 2 LO NN interactions, given in operator form 
times local potentials V(r) in coordinate space. The re- 
production of the NN phase shifts is very good compared 
to the momentum-space N 2 LO NN potentials of Ref. [33| . 
Direct application of the local chiral NN interactions in 
AFDMC sets first nonperturbative benchmarks for the 
neutron matter equation of state at nuclear densities. 
Our results show systematic order-by-order convergence 
with theoretical uncertainties and validate perturbative 
calculations for the softer local NN interactions. Future 
AFDMC calculations with local N 2 LO 3N forces will pro- 
vide ab-initio constraints for nuclear density functionals 
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and for dense matter in astrophysics. This work paves the 
way for QMC calculations with systematic chiral EFT in- 
teractions for nuclei, neutron drops, and nuclear matter, 
can test the perturbativeness of different orders in chiral 
EFT, and allows for matching to lattice QCD results [45 1 
(for example, for few- neutron systems in a box) by vary- 
ing the pion mass in chiral EFT. 
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